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are those who love the entire cosmos... for to them will he shown 
the unbroken web of life and the meaning of infinity. ' 

Alice Walker, The Temple of My Familiar 
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ABSTRACT 



Observations of the microwave sky using the Python telescope in its fifth 
season of operation at the Amundsen- Scott South Pole Station in Antarctica 
are presented. The system consists of a 0.75 m off-axis telescope instrumented 
with a HEMT amplifier-based radiometer having continuum sensitivity from 
37-45 GHz in two frequency bands. With a 0.91° x 1.02° beam the instrument 
fully sampled 598 deg^ of sky, including fields measured during the previous 
four seasons of Python observations. Interpreting the observed fluctuations 
as anisotropy in the cosmic microwave background, we place constraints on 
the angular power spectrum of fluctuations in eight multipole bands up to 
I ~ 260. The observed spectrum is consistent with both the COBE experi- 
ment and previous Python results. Total-power Wiener-flltered maps of the 
CMB are also presented. There is no significant contamination from known 
foregrounds. The results show a discernible rise in the angular power spec- 
trum from large (/ ~ 40) to small (/ ~ 200) angular scales. 
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CHAPTER 1 
INTRODUCTION 

'We know of an ancient radiation that haunts dismembered constellations, a 
faintly glimmering radio station.' 

Cake, Frank Sinatra 

The Cosmic Microwave Background Radiation (CMB), which formed 
when the universe was ~ 300,000 years old (redshift z ~ 1100), is the sig- 
nature of a young, hot, dense universe. Along with big bang nucleosynthesis 
and the Hubble expansion, the CMB is one of the key observational pillars 
that support the big bang model of the universe. 

Since its discovery (Penzias and Wilson 1965) and interpretation (Dicke 
et al. 1965), the CMB has been one of the most powerful tools for discrim- 
inating among cosmological models. The COBE satellite opened up a new 
era in CMB measurement, by detecting anisotropy on the largest scales with 
the DMR instrument (Smoot et al. 1992) and by estabhshing the blackbody 
nature (T=2.728 ± 0.004 K today) of the CMB spectrum to high preci- 
sion with the FIRAS instrument (Mather et al. 1990, Fixsen et al. 1996). 
Measurement of anisotropy in the CMB directly probes conditions of the 
early universe. The anisotropics in the CMB are the seeds of the large scale 
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structures (galaxies, clusters of galaxies, and super-clusters) that we observe 
today. Observations of the angular power spectrum of CMB temperature 
fluctuations can be used to test theories of structure formation and constrain 
cosmological parameters. 

There are three classes of measurements that could be made of the CMB: 
spectral, temperature anisotropy, and polarization anisotropy. This thesis 
describes a specific CMB temperature anisotropy experiment. Python V, with 
an emphasis on data reduction and analysis techniques. It is an expansion 
of Coble et al. (1999). 

1.1 Theory 

The theory of how CMB anisotropics grew via gravitational instability 
from primordial fluctuations is now well understood and has been discussed 
at length by other authors (see, for example, Hu et al. 1997). Before redshift 
z ~ 1100, photons and baryons were tightly coupled via Compton scatter- 
ing, in a photon-baryon fluid. Acoustic oscillations in the fluid were created 
by the resistance of the photon pressure against gravitational compression. 
As the universe cooled, neutral hydrogen formed and radiation decoupled 
from matter and began free-streaming. The epoch of decoupling occurred at 
z ~ 1100 for most cosmologies. Regions of compression and rarefaction in 
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the photon-baryon fluid at decoupling correspond to hot and cold spots that 
we observe in the CMB. 

CMB angular power is usually expressed in terms of C[s. If the sky tem- 
perature is expanded as T{9,(f)) = TiimaimYimi^, 4>)^ then the multipole mo- 
ments, Ci, are defined as C/ =< |a/mP >• Large / corresponds to small angu- 
lar scales. Measuring the angular power spectrum of temperature anisotropics 
in the CMB is a powerful tool for determining many cosmological parame- 
ters, including the expansion rate of the universe {Ho), the total density of 
the universe {Q), the amount of baryonic matter {Qf,) in the universe, the 
amount and nature of dark matter, the spectral index of primordial fluctu- 
ations (n), the cosmological constant (A), and the ionization history of the 
universe. Figure |I| shows the dependence of the angular power spectrum on 
various cosmological parameters. 

Traditional cosmological tests, such as galaxy and cluster distributions, 
cannot probe these parameters directly, because they only determine the dis- 
tribution of visible mass. They must assume an underlying bias, the unknown 
relation between mass and light, whereas CMB anisotropy corresponds di- 
rectly to density fluctuations (White et al. 1994). CMB measurements are 
also independent of the extra-galactic distance ladder and probe parameters 
on a global scale. 
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As one example of the dependence of the CMB angular power spectrum 
on a cosmological parameter, the angular scale of the first acoustic peak 
(Figure |I[) is a sensitive indicator of the total density, Q. The horizon scale 
at the epoch of decoupling, which corresponds to the scale of the first acoustic 
peak, is given by 

AO ^ 0.87°Q^/\ (1.1) 

As Q decreases, the peak occurs at smaller angular scales. For a flat universe 
{Q=l), the corresponding angular scale is ~ 1°. 

1.2 Temperature Anisotropy Experiments 

As with most astronomical measurements, CMB observations are usu- 
ally conducted from locations inhospitable to humans, namely places that 
are high, dry, and cold, in order to reduce the interference from the atmo- 
sphere. CMB experiments have been and will be conducted from satellites, 
such as COBE, MAP (see the MAP website), and Planck (Planck website). 
The advantage to satellite missions is that they can map the whole sky with 
high precision and without interference from the atmosphere. The disadvan- 
tages are that they often have a long development period and are costly. A 
number of balloon-borne experiments have been conducted, such as ARGO 
(de Bernardis et al. 1994, Masi et al. 1996), BAM (Tucker et al. 1996), 
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Figure 1 Model CMB angular power spectra for various cosmological parame- 
ters. Large / corresponds to small angular scales. Shown are a standard CDM 
model (solid line, density = 1.0, Hubble constant Hq = 50 km/s/Mpc, 
baryon density Qi, — 0.05, and cosmological constant A = 0.0), a model 
with high baryon content (dot-dashed line, Q — 1.0, Hq — 50 km/s/Mpc, 
Qb = 0.20, and A = 0.0), a model with a non-zero cosmological constant 
(dotted line, Vl = 1.0, Ho = 50 km/s/Mpc, = 0.02, and A = 0.70), and 
an open model (dashed line, = 0.26, Hq — 50 km/s/Mpc, = 0.06, 
and A = 0.0). The models were computed using CMBFAST (Seljak and 
Zaldarriaga 1996). 

FIRS (Ganga et al. 1994), MAX (Lim et al. 1996, Tanaka et al. 1996), 
MSAM (Cheng et al. 1997), and QMAP (Devhn et al. 1998, Herbig et 
al. 1998). These balloon flights typically observe the sky for one night per 
year. The next generation of balloon-borne experiments, long-duration bal- 
looning (LDB) experiments, are currently under way. The BOOMERANG 
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(BOOMERANG website) experiment circumnavigated Antarctica, observing 
for approximately 10 days in December 1998. The TopHat (TopHat website) 
experiment is scheduled for an LDB flight in December 1999. A plethora of 
experiments have detected or constrained CMB anisotropy from the ground. 
The South Pole is an especially attractive site for CMB measurement be- 
cause of its stable atmosphere and high altitude. The Python (Dragovan et 
al. 1994, Ruhl et al. 1995, Piatt et al. 1997, Kovac et al. 1999), Viper 
(Peterson et al. 1998), and ACME (Gaier et al. 1992, Schuster et al. 1993, 
Gundersen et al. 1995) telescopes have made observations from the South 
Pole. Other ground-based experiments include CAT (Scott et al. 1996), lAB 
(Piccirillo and Cahsse 1993), lAC (Femenia et al. 1997), MAT (Torbet et al. 
1999), OVRO NCP (Readhead et al. 1989), OVRO RING5M (Leitch et al. 
1998), Saskatoon (SK) (Netterfield et al. 1996), SuZIE (Church et al. 1997) 
and Tenerife (Hancock et al. 1997). 

Since anisotropics are temperature differences, CMB experiments must 
difference (chop) in some way. The chopping, whether internal or external, 
must be quick, and come back to the same place on the sky faster than 
the 1/f noise in the electronics and atmosphere. Often this is accomplished 
with a chopping secondary mirror (the chopper) and/or an internal Dicke 
switch. Some simple chopping strategies include a single difference (2-point 
chop), for which ATj = Tj+i— Tj, a double difference (3-point chop), for which 
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ATi ^Ti- - ir^+i, and a 4-pomt chop, for which AT^ = -1^^+ fTl+i - 
f + \Ti+z- Other experiments perform a continuous (fast) scan with the 
chopper and modulate the data in software, scan slowly with the telescope 
and read out the detectors quickly, or spin rapidly. Chopping strategies are 
sometimes called lockins, harmonics, modulations, or demodulations. 

There are 2 basic detector technologies used for CMB experiments: bolome- 
ters and HEMTs. SIS detectors are also sometimes used. HEMTs have the 
advantage that they are robust and easy to use. They are inherently polar- 
ization sensitive, coherent detectors that typically operate at frequencies ^ 
90 GHz. They can operate at room temperature, but become more sensi- 
tive as they are cooled down to ~ 10 K. Bolometers are more sensitive than 
HEMTs, but they are less robust and must be cooled to at least 0.25 K, so 
much work must go into the cryogenics. Bolometers typically operate at fre- 
quencies ^ 90 GHz and are often used on balloon-borne experiments, where 
their sensitivity is better than the atmospheric noise. It is advantageous to 
make measurements at several frequencies in order to spectrally distinguish 
possible contamination by galactic foregrounds (e.g. dust, synchrotron). 

CMB anisotropics can also be measured using interferometers, such as 
CAT, which has detected anisotropy, DASI (White et al. 1997), which will 
observe at large to medium angular scales from the South Pole, and CBI 
(CBI website), which will observe at small angular scales from Chile. Inter- 
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ferometers are advantageous because they have low systematics. 

In choosing an observing strategy (how much time to spend on which 
region of the sky and in what order), the theoretically optimal signal-to-noise 
ratio S/N=l, is too optimistic. In reality, experiments should have S/N>1, 
in order to perform consistency checks on the data. Interlaced observing 
strategies are desirable because they reduce striping in the maps obtained, 
but that is not always possible given observing sites and conditions. 

Results from the COBE satellite (Smoot et al. 1992) tightly constrain the 
angular power spectrum at the largest angular scales. Several experiments 
have measured the angular power spectrum at degree angular scales (e.g., 
SP, MAX, Python I-IV, MSAM, SK, QMAP). Figure | is a compilation of 
measurements of the angular power spectrum previous to the publication of 
the PyV measurements. Collectively the data show a rise in power towards 
smaller angular scales; individually no experiment covers a wide range of 
angular scales from COBE scales to degree scales, and most cover only small 
regions of the sky. 

The dataset from the fifth Python observing season (hereafter PyV) has 
sufficient sky coverage to probe the smallest scales to which COBE was sen- 
sitive, while having a small enough beam to detect the rise in angular power 
at degree angular scales. 



9 



100 



80 



60 



40 



20 




+ COBE 
m FIRS 

Tenerife 

SP 

BAM 

□ ARGO 
K MAX 

Python 
lAC 
lAB 

m QMAP 

O SK 

CAT 

OVRO RING5M 

□ SuZIE 

□ MSAM 
OVRO NCP 



10 



100 

Multipole 1 



1000 



Figure 2 CMB anisotropy measurements as of Feb 1999. The points are 
COBE (Hinshaw et al. 1996), FIRS (Ganga et al. 1994), Tenerife (Hancock 
et al. 1997), SP (Gundersen et al. 1995), BAM (Tucker et al. 1996), ARGO 
(de Bernardis et al. 1994, Masi et al. 1996), MAX (Tanaka et al. 1996), 
Python I-III (Piatt et al. 1997), TAB (Piccirillo and Calisse 1993), lAC 
(Femenia et al. 1997), QMAP (de Oliveira-Costa et al. 1998), SK (Netterfield 
et al. 1996), CAT (Scott et al. 1996), OVRO RING5M (Leitch et al. 1998), 
SuZIE (Church et al. 1997), MSAM (Wilson et al. 1999), and OVRO NCP 
(Readhead et al. 1989). 



CHAPTER 2 
THE PYTHON EXPERIMENT 



'Python is a large non- venomous Old World snake that coils around and suf- 
focates its prey.' 

Oxford American Dictionary 

2.1 Introduction to the Python Experiment 

In its first four seasons the Python experiment detected significant anisotropy 
in the CMB (Dragovan et al. 1994 (Pyl), Ruhl et al. 1995 (Pyll), Piatt et 
al. 1997 (Pylll), Kovac et al. 1999 (PylV)). Observations from the first 
three seasons were made at 90 GHz with a bolometer system and a 4-point 
chop scan strategy, yielding CMB detections at angular scales of / ~ 90 and 
I ~ 170. During the PylV season measurements were made using the same 
scan strategy with a HEMT amplifier-based radiometer, confirming Pyl-III 
detections in a 37-45 GHz frequency band. 

Observations were made from November 1996 through February 1997 in 
the fifth Python observing season. In order to increase the range of observed 
angular scales, a smoothly scanning sampling scheme was implemented. As 
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Season 


dates 




FWHM 


chop 


throw 


^eff 


fields 


Pyl 


12/92 


90.3 


0.82° 


4-pt 


8.25° 


87 


15 


Pyll 


12/93 


90.3 


0.82° 


4-pt 


8.25° 


87 


31 


Pylll-L 


10/94-12/94 


90.3 


0.82° 


4-pt 


8.25° 


87 


110 


Pylll-S 


10/94-12/94 


90.3 


0.82° 


4-pt 


2.75° 


170 


140 


PylV 


2/96 


40.1 


1.06°xl.l2° 


4-pt 


8.25° 


85 


14 


PyV 


11/96-2/97 


40.3 


0.9rxl.02° 


triangle 


11.0° 


40-260 


598 deg2 



Table 1 Python observing parameters (adapted from Kovac et al. 1999). Uc is 
the central frequency of the observations, FWHM is the FWHM of the beam, 
chop is the chopping strategy, throw is the extent of the chopper sweep, leff 
is the multipole (angular scale) at which observations were made, and fields 
is the number of fields observed or the sky coverage. 

a result, PyV is sensitive to the CMB angular power spectrum from / ~ 40 
to / ~ 260. 

Table ^ summarizes the observing parameters from all five Python sea- 
sons. 



2.2 Instrument 



The PyV measurements were made using the same receiver as the PylV 
system as described in Alvarez (1996) and Kovac et al. (1999). A schematic 
of the receiver is shown in Figure ^. The receiver consists of two focal-plane 
feeds, each with a single 37-45 GHz HEMT amplifier. A diplexer splits each 
signal at ~ 41 GHz before detection, giving four data channels. The analysis 
reported here eventually combines signals from all four channels, resulting in 
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Figure 3 Schematic diagram of the receiver (Alvarez 1996). The signal path 
is shown for only one of the two feeds. 

a thermal radiation centroid i>c = 40.3 GHz and effective passband Ai/ = 6.5 
GHz for the PyV dataset. The bandpasses for the four data channels are 
shown in Figure ^. 

The receiver is mounted on a 0.75 m diameter off-axis parabolic telescope 
(Dragovan et al. 1994), which is surrounded by a large ground shield to 
block stray radiation from the ground and Sun. The beams corresponding 
to the two feeds observe the same elevation and are separated by 2.80° on 



13 
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Figure 4 Bandpasses for the 4 data channels. Channel 3 is a dark channel. 
Channels 1 and 2 observe the same points on the sky as do channels 4 and 5. 
The bandpasses are similar in frequency range, so that the reduced data from 
channels which observe the same points on the sky are eventually co-added 
(Chapter 3). 

the sky. These beams are scanned horizontally across the sky by a large 
rotating vertical flat mirror, the chopper, at 5.1 Hz. The new scan strategy 
motivated two changes from the instrument configuration described in Kovac 
et al. (1999): the frequency response of the data system was extended by 
switching to 100 Hz-rolloff antialiasing Bessel filters, and the data recording 
rate was correspondingly increased, to 652.8 samples/sec for each channel. 
A schematic of the telescope is shown in Figure ^ Pictures of the telescope 
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with and without the ground shield are shown in Figures |^ and |^ respectively. 




Figure 5 Schematic diagram of the Python telescope (Dragovan et al. 1994). 
The telescope is surrounded by a large ground/Sun shield. The chopping flat 
(chopper) is oriented vertically and scans in azimuth. 

2.3 Calibration 

As in previous Python seasons, the primary DC calibration of the detec- 
tors was derived using liquid nitrogen, liquid oxygen, and ambient temper- 
ature thermal loads external to the receiver (Dragovan et al. 1994, Ruhl et 
al. 1995). Load calibrations were performed approximately once per day, 
and gains were found to be consistent over the entire season to within ± 




Figure 6 Picture of the Python telescope without the ground shield or dewar. 
For CMB observations the telescope is oriented such that the chopper is 
vertical. 
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Figure 7 A view of the Python telescope and surrounding labs. 

2%, with no discernible trends. Gain compression, which was a source of 
systematic uncertainty in the calibration of Pyl-Pylll, is measured to be 
negligible for the Python HEMT receiver. Systematic uncertainty in the DC 
load calibration is estimated to be ± 10%. 

Several efficiencies must be estimated to relate the load calibrations to 
celestial response in the main beam, which account for power losses in the 
atmosphere, in the sidelobes, and in the telescope, and they are calculated 
using data from skydips and from various beam measurements. In order to 
account for the net effect of atmospheric attenuation and beam efficiency, the 
data must be multiplied by a factor of 1.10. The resulting systematic calibra- 
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tion uncertainty of ^4%° is asymmetric, due to the fact that the individual 
losses are small positive numbers and hence the errors in their estimation 
follow skewed distributions (Kovac et al. 1999). 

The dynamic response of the system was calculated from laboratory mea- 
surements of the transfer functions for the AC coupling and antialiasing fil- 
ters in the data system, and confirmed on the telescope by comparison of 
observations made of the moon using normal and slow chopper speeds. The 
response speed of the detectors is not a concern for this calibration or for 
its uncertainty. An appropriate response correction factor is applied to each 
modulation of the data. The uncertainty on these factors is small, and is 
dominated by a ± 5% systematic uncertainty on their common normaliza- 
tion. 

The overall uncertainty in the calibration of this dataset is estimated to be 
Antenna temperature has been converted to units of STqub through- 
out. The calibrations, efficiencies and their uncertainties are summarized in 
Table |. 

2.4 Observations 

Observations for the PyV season were taken from November 1996 through 
February 1997. 
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Calibration 



Uncertainty (%) 



DC load systematic 
DC load statistical 
atmospheric attenuation 
plus beam efficiency 
dynamic response 
TOTAL 



± 10 
± 2 



+10, -4 
± 5 

+15, -12 



Table 2 Calibrations, efficiencies and corresponding uncertainties. 
2.4-1 Observing Regions 

Two regions of sky were observed: the PyV main field, a 7.5° x 67.7° 
region of sky centered at a = 23.18'^, 5 = —48.58° (J2000) which includes 
fields measured during the previous four seasons of Python observations and 
a 3.0° X 30.0° region of sky centered at a = 3.00'^, 6 = -62.01° (J2000), which 
encompasses the region observed with the ACME telescope (Gundersen et 
al. 1995). The total sky coverage for the PyV regions is 598 deg^, greater 
sky coverage than previous degree-scale CMB experiments. Figure ^ shows 
the IRAS 100 nm view of the galaxy with the Python observing regions 
superimposed. 

2.4-2 Observing Strategy 

Both PyV regions are fully sampled with a grid spacing of 0.92° in eleva- 
tion and 2.5° in right ascension, corresponding to a distance of 1.6° on the 



19 




Figure 8 IRAS 100 fim view of the galaxy with the Python observing re- 
gions superimposed. The 2 large white boxes indicate the PyV regions, the 
medium- sized orange box Pylll, and the small red box PylV. This demon- 
strates the large sky coverage achievable from the South Pole in one of the 
most clean regions of sky. 

sky at a declination of —50°. A grid of the center positions of the PyV fields 
is shown in Figure ^. The telescope is positioned on one of the fields and the 
chopper smoothly scans the beams in azimuth in a nearly triangular wave 



pattern (Figure [T^). The chopper throw is 17° in azimuth, corresponding to 
11° on the sky at a declination of —50°. The telescope tracks on each field; 
it does not drift scan. 



There are 128 data samples for each detector channel in a complete chop- 
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Figure 9 Center positions of the Python V fields. The fields are 2.5° apart 
in right ascension, which corresponds to 1.6° on the sky at a declination of 
—50°. The fields are 0.92° apart in declination. 

per cycle, and 164 chopper cycles (~ 32 s) of data (one stare) are taken of 
a given field before the telescope is positioned on the next field in the set. 
One data file consists of 164 chopper cycles for each field in the set. For a 
set which contains 13 fields, one data file is ~ 32 s x 13 ~ 6.9 min long. 

Approximately 13 hours of good data (100 files) are taken of a set of 
fields before the telescope moves on to the next set of fields. A total of 309 
fields are observed in 31 sets of 5-17 fields. Some fields are observed in more 



than one set, yielding 345 effective fields. Figure |Tl] labels the observing sets 
with respect to the grid shown in Figure ^ and Table ^ lists the number of 
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Figure 10 Motion of the chopper. The chopper smoothly scans in azimuth 
in a nearly triangular wave pattern at 5.1 Hz. The waveform is not perfectly 
triangular because the chopper slows down when turning around and because 
the LVDT is non-linear. 

fields in each set, the number of files of good data taken in each set, and the 
declination at which each set was taken. 



2.4.3 Pointing 

The combined absolute and relative pointing uncertainty is estimated to 
be 0.15°, as determined by measurements of the moon and the Carinae nebula 
{a = 10.73/i, (5 = —59.65°). The beams corresponding to the 2 feeds observe 



the same elevation and are separated by 2.80° on the sky (Figure |T^). 
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Figure 11 Observing sets. Each set is observed for approximately a day (leav- 
ing ~ 13 hours of good data) before moving on to the next set. Neighboring 
sets in the main PyV region overlap by 1 field. Circles represent the fields 
and arrows point to the end fields in each set. In addition to the sets shown, 
sets ib, jb, kb, and lb were observed with a scan pattern of 5 fields per file. 
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Figure 12 Pointing of the two feeds on the sky in the main PyV region. 
The beams corresponding to the two feeds observe the same elevation and 
are separated by 2.80° on the sky. The squares represent the positions of 
channels 1 and 2 and the circles represent the positions of channels 4 and 5. 
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2.44 Beam 

The beam is determined from scans of the Carinae nebula and the Moon. 
Scans of the Carinae nebula were obtained from the Viper telescope (Griffin 
and Peterson 1998), which has the same receiver system as PyV and a 0.25° 
beam. The Viper scans were convolved with asymmetric Gaussian beams 
of widths ax and ay. The widths were varied as free parameters and best 
fits to the observed PyV beams were obtained. This convolution procedure 
works better than a deconvolution procedure because of the noise in the 
measurements. The results are not dependent on size of the Carina nebula 
assumed in the fits. The PyV beam is well approximated by an asymmetric 
Gaussian of FWHM 0.911^? x 1-021^? degrees (az x el). 
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Set 


-^fields 


Affiles 


5 


sa 


17 


41 


-62.99 


sb 


17 


103 


-62.01 


sc 


17 


33 


-61.01 


ga 


7 


302 


-51.80 


gb 


9 


184 


-51.80 


gc 


13 


95 


-51.80 


ha 


13 


100 


-50.88 


hb 


9 


84 


-50.88 


he 


13 


104 


-50.88 


ia 


13 


45 


-49.97 


ib 


9 


130 


-49.97 


ib 


5 


52 


-49.97 


ic 


13 


119 


-49.97 


ja 


13 


120 


-49.03 


jb 


9 


269 


-49.03 


jb 


5 


250 


-49.03 


ic 

J 


13 


109 


-49.03 


ka 


13 


81 


-48.14 


kb 


9 


186 


-48.14 


kb 


5 


295 


-48.14 


kc 


13 


59 


-48.14 


la 


13 


191 


-47.20 


lb 


9 


298 


-47.20 


lb 


5 


343 


-47.20 


Ic 


13 


122 


-47.20 


ma 


13 


113 


-46.28 


mb 


9 


198 


-46.28 


me 


13 


124 


-46.28 


na 


13 


153 


-45.36 


nb 


9 


50 


-45.36 


nc 


13 


107 


-45.36 



Table 3 Observing sets. 



CHAPTER 3 
DATA REDUCTION 

3.1 Data Cuts 

A total 719.52 hours of data were taken during the PyV season. After 
cutting 45.6% of the data for file errors, tracking errors and weather, 391.32 
hours of data remained for use in the CMB analysis. Table § summarizes the 
data cuts. 

A file is cut: 

• for file errors if any of the data, pointing, or time matrices are empty, 
if the time matrices from the two data computers did not agree, or if 
the number of stares in the data, pointing and time matrices did not 
agree. 

• for a relative pointing error if the stare that is furthest from where 
it should be relative to the scan center is more than 3' from where it 
should be. 

• for an absolute pointing error if its scan center deviates by more than 
6' from the scan center of all other files in the set. 
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Cut 



hours percent 



internal file errors 
relative pointing 
absolute pointing 
chopper error 



10.42 1.45 
49.39 6.86 

73.43 10.21 
7.47 1.04 
24.48 3.40 
163.01 22.65 
328.20 45.61 



low levels 
high levels 



TOTALCUTS 



Table 4 Data cuts. 

• if the mean chopper waveform in any stare deviates by more than 1% 
from the benchmark waveform made from the whole data set. 

• for anomalously low noise levels, which could be caused by the power 
failing and the dewar heating up or the gain being set to the wrong 



• for anomalously high noise levels, which are caused by weather. There 
is a clear distinction between good and bad weather. 



3.2 Modulations 



The data are modulated using 



level. 



1 



m 



= 1 



cosim-KO/Oc) X < 



(3.1) 



H{e) m = 2...8 
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where m is the modulation number, 6 is the chopper angle, 6c is the extent 
of the chopper throw and H{6) = 0.5(1 — cos(2tt6/6c)) is a Hann win- 
dow (Figure 0). The m = 2 . . . 8 modulations are apodized with the Hann 
window in order to reduce the ringing of the window functions in Z-space. 
Data taken during the right- and left-going portions of the chopper cycle are 
modulated separately, to allow for cross-checks of the data. Sine modulations 
are not used in the analysis because they are anti-symmetric and are thus 
sensitive to gradients on the sky. 
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Figure 13 The data are modulated with cosines which have been apodized 
by a Hann window. Three of the modulations are shown. 



The modulations in equation 



3.1 



and Figure ^ are cosines or apodized 
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cosines in space, not in time. Since the chopper waveform (Figure ^) is 
non-hnear, the timestream modulations are not strictly cosines or apodized 
cosines in time (or chopper sample). Rather, they are created to produce the 



spatial response in equation and Figure |13[ The timestream modulations 



as a function of chopper sample are shown in Figure The timestream 
modulations, Mm(s) are normalized such that J^s^ = (so that a uniform 
temperature does not contribute) and J2s \ M\ = 2 where s is chopper sample, 
following the single difference convention. 
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Figure 14 Timestream modulations. Since the chopper waveform is non- 
linear, the timestream modulations are created so that the spatial response 
will be cosines or apodized cosines. These 3 timestream modulations corre- 



spond to the 3 spatial modulations shown in Figure [13. 
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3.3 Instrument Transfer Function 



The instrument transfer function produces two main effects on the data: 
a roll-off with frequency and a phase shift between data taken during the 
right and left-going portions of the chopper cycle (hereafter right and left 



data). Figure |T5| shows the frequency dependence of the transfer function. 
Since each modulation is well-isolated in frequency space, the effect of the 



1.2 
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Figure 15 Frequency dependence of the instrument transfer function. 

frequency roll-off of the transfer function can be accounted for with 1 number 
for each modulation (Table §). The modulated data are divided by T to 
remove the frequency dependence of the transfer function. 



Figure |16| shows the phase effect of the transfer function on right and 
left unmodulated moon data and Figure O shows modulated moon data, in 
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mode 


/ (Hz) 


T 


1 


13 


1.0 


2 


22 


.99 


3 


33 


.98 


4 


43 


.95 


5 


55 


.92 


6 


66 


.88 


7 


77 


.83 


8 


88 


.79 



Table 5 Instrument transfer function vs. mode. Frequency / is the peak 
frequency of each modulation. 

which right and left data are properly phased. If the right and left data are 
not properly phased, the signal will be washed out. 
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Figure 16 Unmodulated moon data showing effect of the transfer function 
on the right-left phase. Data from one stare and channel is shown as an 
example. 
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Figure 17 Checking the phase of right and left data on the moon. Data from 
channel 1, modulation 3 is shown as an example. If the right and left data 
are not properly phased, the signal will be washed out. 

3.4 Chopper Synchronous Offset 



Data in a given file, field, channel, and modulation are co-added over all 
chopper cycles. A chopper synchronous offset, caused by differential spillover 
past the chopper, is removed from each data file by subtracting the average 
of all of the fields in a file. This is not just a DC offset; there is an offset 
removed for each modulation (or equivalently, for each sample of the chopper 



waveform) and channel. Figure shows the chopper offset as a function of 



time for one of the channels and modulations as an example. Figure shows 
a histogram of the chopper offset for that same channel and modulation. The 



large scatter in both figures is due to noise. 
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Figure 18 Chopper offset vs. time for channel 1, modulation 2, right-going 
data. Each point corresponds to the offset removed from one file. 

To test that the chopper offset is stable over a file, we take the difference 
in chopper offset between two adjacent files in time. A plot of the chopper 
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Figure 19 Histogram of chopper offset for channel 1, modulation 2, right- 
going data. 
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Figure 20 Chopper offset difference vs. time for channel 1, modulation 2, 
right-going data. Each point corresponds to the offset difference between 
adjacent files in time. 

offset difference is shown in Figure |2^ and a histogram is shown in Figure 



21. Again, there is a large scatter due to noise. 
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Figure 21 Histogram of chopper offset difference for channel 1, modulation 
2, right-going data. 
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In order to compare the size of the chopper offset and the chopper off- 
set difference with the scatter due to noise reduced, we compute the mean 
chopper offset and mean chopper offset difference for each channel and mod- 



ulation over the whole data set. As illustrated in Figures and |23|, the 
chopper offset difference is significantly smaller (by at least 3 orders of mag- 
nitude) than the chopper offset. Further, when the differencing procedure is 
repeated as a function of file lag, the results indicate that the chopper offset 
is stable over ~ 5 files. Thus our assumption that the chopper offset is stable 
over 1 file is a correct and even a conservative assumption. 
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Figure 22 Mean chopper offset over the data set for all channels, right and 
left-going data vs. modulation. A different chopper offset is subtracted from 
each file. 
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Figure 23 Mean chopper offset difference over the data set for all channels, 
right and left-going data vs. modulation. The chopper offset difference is 
significantly smaller that the chopper offset (Figure ^2[), indicating that the 
chopper offset is stable over a file. The vertical scale in this figure covers a 
smaller range than that of Figure E2. 



3.5 Ground Shield Offset 



When the modulated data are binned in azimuth, a periodic signal due 
to the 12 panels of the ground shield is evident, especially on larger angular 
scales. The signal of period 30° is fit for an amplitude and is subtracted. 
Removal of the ground shield offset has less than 4% effect on the final angular 
power spectrum in all modulations, because when the data are binned in RA, 
the effect averages out. 

The ground shield offset is computed from all fields and files in an observ- 
ing set. A different ground shield offset is fit and removed for each channel. 
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modulation, and set. Figure ^ shows the azimuthally binned data for one 
set, the fit to the ground shield synchronous signal, and the same data with 
the ground shield offset removed. The ground shield offset is larger for data 
taken when the telescope is pointed closer to the ground shield and for mod- 



ulations which probe larger angular scales, as illustrated in Figure |25 
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Figure 24 Top- Azimuthally binned data for channel 1, left-going chopper 
data, modulation 1, set jb9. There are 2421 data points (269 files x 9 fields) 
plotted. The periodic signal caused by the 12 panels of the ground shield is 
evident and the best fit curve is plotted as a solid line. Bottom- the same 
data with the ground shield offset subtracted. 
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Figure 25 Amplitude of the ground shield offset as a function of declination 
and modulation for channel 1, left-going chopper cycles, file sets sb, jb, and 
nb. The offset is larger for data taken at a lower elevation, i.e. when the 
telescope is pointed closer to the ground shield. 

3.6 Co- adding 

After the data have been modulated and offsets removed, the right and 
left-going data, which have been properly phased, are co-added, as are data 
from channels which observe the same points on the sky. The bandpasses 
are similar enough in frequency range that the data can be co-added without 
loss of information. Hereafter channel 12 will refer to the co-added chan- 
nel 1 and 2 data and channel 45 will refer to the co-added channel 4 and 5 
data. Since the two feeds observe different points on the sky they cannot be 
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co-added; the theoretical and noise covariances between them are included 
in the likelihood analysis of the angular power spectrum. 

3.7 60 Hz 

The South Pole Station power nominally operates at a frequency of 60 Hz, 
so there is possible contamination at 60 Hz and its harmonics. However, the 
chopper runs at 5.1 Hz, which is incommensurate with 60 Hz. To verify that 
a 60 Hz signal would not contaminate the data, we created a 60 Hz signal 
for one stare of data, modulated it, averaged over 164 cycles and found that 
the resulting signal is negligible for all modulations. 



CHAPTER 4 
DATA ANALYSIS TECHNIQUES 

4.1 Bayesian Likelihood Analysis 

In order to estimate parameters from data, a Bayesian likelihood analysis 
(see for example Readhead et al. 1989) is used. What we want to report 
from any experiment is the probability of the theory being true given the 
data, P {theory \ data) . What we can actually measure is the likelihood, C — 
P{data\theory) . Applying Bayes' theorem, 

P {theory \data) oc P {data\theory)P {theory) . (4.1) 

The prior probability, P {theory), is conservatively taken to be uniform (un- 
informative) . This technique assumes the data are drawn from a Gaussian 
distribution on the sky. 

The likelihood {C) can be written: 

C = (27r)-^/2det(C)-'/'exp(-xV2) (4.2) 

where = SC-^d, = 0^ + 0^ + C^, cTis the data, is the theory 
covariance matrix, C'^ is the noise covariance matrix, and C*-^ is the con- 
straint matrix for offset subtraction. for PyV is discussed in section 4.4 

and C'^ for PyV is discussed in sections 5.2 and 6.1. 
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Figure 26 Theory covariance matrix, C"^, for PyV modulation 2. C"^(l, 1) 
is in the bottom left. The matrix is 690 x 690 for the 345 fields and 2 
channels in the reduced, co-added data. Points 1:345 correspond to channel 
45 and points 346:690 correspond to channel 12. There are significant theory 
correlations between channels because they are separated by 2.80° on the 
sky. 

The theory covariance matrix, C^, depends on the model CMB angu- 



lar power spectrum Ci and the experimental window functions Wuj. The 
elements of C'^ are given by 



4:71 



iVVlij- 



(4.3) 



A sample theory covariance matrix is shown in Figure |2B. 
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4.2 Window Functions 

The window functions are a measure of experimental sensitivity as a func- 
tion of angular scale /. They are generated from the experimental beam map, 
modulations, and observing strategy. 

4.2.1 A Method For Calculating Window Functions 

The following is a method for calculating the window functions for a 
general beam map and observing strategy in the fiat sky approximation. 
Section 4.2.3 discusses the importance of using the exact beam map and 
observing strategy for certain experiments. The flat sky approximation is 
valid for many current and past CMB experiments. Another discussion of 



window functions can be found in White and Srednicki (1995). Figure |27 
illustrates the geometry of the beam map and observing strategy used in the 
derivation below. 

• The (modulated) beam map is B, which is 2-dimensional. Bm{x — Xi) 
is the beam map for modulation m corresponding to an observation at 
pixel i. 

• 9i is the twist of the beam map for pixel i with respect to the Si-axis 
in Figure p7. 
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Figure 27 Position of beam maps for computing Wuj. si and S2 are flat space 
coordinates. 



The rotation matrix acting on [x — Xi) is Ri 



/ \ 

cos{9i) —sin{9i) 



^ sin{9i) cos{9i) ^ 

• The measured anisotropy at a point is a convolution of the true anisotropy 
with the beam map: ATj = / Bm^R-iix — Xi))AT{x)dx. 

• The underlying temperature field AT{x) can be expanded in spherical 
harmonics: AT(x) = Ili^aimYim{9 , (p) , and the multipole moments, C;, 
are Q =< \aim\^ >. 



The elements of the theory covariance matrix are Cfj — (ATiATj) 
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Using the above definitions, the elements of the theory covariance matrix 



are: 



J J BrniRiix- x^))B.m'(Riix - Xj)){AT{x)AT{x))dxdx. (4.4) 



Assuming rotational symmetry and averaging over the whole sky, 



(AT(f)Ar(f' )) = E ^^^QPiicosiq)), (4.5) 



where Pi are the Legendre polynomials associated with the spherical har- 
monic expansion, q = \q\, and 

q = X — x' . (4.6) 

Defining dummy variables, 

V = Rifx — x.?\ 

(4.7) 

f = Rj(f'-fj), 

we have 



y = 




y' = 




X = 


Ri~\y)+Xi 


-V 

X = 


Rf\y')+x, 


dx = 


Rr\dy) 


dx' = 


Rf'idy') 



(4.8) 



Substituting these into equation we have 



Cl,=l lBUy)Brni^)T.^^^CiPiicosiq))Rr\dmr\dy') (4.9) 
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Defining the window functions Wuj by 



(4.10) 



and combining equations [4.9| and [4.10| , we get an expression for the window 
functions: 



Wu, = J J BUy)Bm'iy')Piicosq)Rr\dy}Ri 



(4.11) 



This is a four-dimensional integral, which is computationally intensive. 
In order to reduce the dimensionality of the integral to be able calculate the 
window functions in a reasonable amount of time on a workstation, we make 



the fiat sky approxmation that for large 1, 



Pi{cosq) = Jo{lq). 



(4.12) 



The expression for the window functions then becomes 



Wh, = J J BUy)Bn.'{y')Mlq)Rr\dy)Rf\dy'). (4.13) 



Using the identity 



Joilq) = (l/27r) / d(Pe 



-ilqcostf) 



(4.14) 



and substituting into equation 4.13, we have 



Wuj = (l/27r) / / / BUy)BrrAy')e-''''°''^d(j)Rr\dy}Rf\dy'). (4.15) 



Defining a vector k such that | A;| = / and k-q = Iqcoscp, equation [4.15| becomes 



Wk^j = (l/27r) / / / Bm{y)Bmiy')e-''^-'^d<pRr\dy)Ri~\dg). (4.16) 



45 



From equations |4.6| and |4.8|, we have 



(4.17) 



Substituting equation [4.17| into equation |4.16|, the window functions become 



Wmj = (l/27r) y y y i?^(y)i?^,(y')e-^'=-(^'" fi7)+-.-Rr (2^)-^.)ci0Rr^((iy)Rj-^(ciy'). 

(4.18) 

Since we are integrating over all y and since the determinant of the rotation 
matrix is 1, we can write this as 

Wkij = (l/27r) III 5^(y)5^,(^)e-*^^-(^'"(^^)+^^-^j"(3r)--.)c^0rf^rf^. 

(4.19) 

Rearranging a few terms, 



W,,, = (l/27r) / e-''^-^^-^- 



(4.20) 



Equation 4.20 can be written in terms of the Fourier transformed beam maps, 



yielding the final expression for the window functions: 



Wkij = (l/27r) / e-*^^(^^-^^^)E^(fcRi-i)5;;,(fcRj-i)rf0. (4.21) 



The four-dimensional integral of equation |4.11| has been reduced to a 
one-dimensional integral and a Fourier transform, saving computing time. 
The window functions can be saved to disk for quick computation of C'^ for 
various input Ci spectra using equation [4.10. 



We will refer to the array Wuj as the window functions and to Wi 
as the diagonal window functions. 
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= Wui 



4-2.2 Normalization of the Beam Map 

The modulated beam map is created by convolving the (spatial) modula- 
tion pattern with the beam and is usually normalized by (l/27rcr^), where a 



characterizes the width of a symmetric Gaussian beam, or a = y-^icTxC^y) in 
the case of an asymmetric Gaussian. The normalization of the modulations 
usually follows the convention discussed in 3.2. 

Some experiments do not follow this beam normalization. For example, 
for the MSAM-I experiment [(Cheng et al. 1997)| , which does not have a 



Gaussian beam, the beam map and data are normalized consistently with 
each other, but the data are in units of nKdeg^ rather than nK. 

4-2.3 Approximations 

Depending on the experiment, some approximations can affect the win- 
dow functions. For example, for the MSAM-I experiment, a Gaussian beam 
is a poor approximation. Its impact on the window functions is shown in 



Figure ESI. Ignoring the twist of the MSAM-I experiment is also a poor 



approximation, as shown in Figure 29 
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Figure 28 Diagonal window functions for MSAM-I with a Gaussian beam 
and with the true beam. 

4.2.4 Analytic Expressions 

This section provides some concrete examples of window functions for 
observing strategies and beam maps which can be treated analytically. Ex- 
amples will be given for 3 cases of analytic expresssions: (1) calculation of 
the full window function array in the flat space approximation, (2) calcula- 
tion of the diagonal window functions in the flat space approximation, and 
(3) calculation of the diagonal window functions without the flat space ap- 
proximation. 



MSAM Beam 
Gaussian Beam 
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Figure 29 Effect of twist on tlie window functions. Including twist is impor- 
tant for tlie MSAM-I experiment. 

• Case 1: Two examples of the full window function array given by equa- 
tion |4.21| are given below. 



Example 1.1. Gaussian beam. 

The beam map for a symmetric Gaussian beam of width a is given by 



so the Fourier transform of the beam map is 



(4.22) 



B{k) = e 



(4.23) 
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Using equation [4.21| the window functions are: 



(4.24) 



Jo{l{\Xi - Xj\). 



Example 1.2. Single difference Gaussian beam. 

The beam map for a symmetric Gaussian beam modulated with a single 



difference chop is given by 



Bio] 



27ra2' 



,-|a-ao|V{2T") _ p-|«+ao|V(2cr^) 



(4.25) 



and its Fourier transform is 



(4.26) 



Using equation [4.21| , the window functions are 
Wkij 



1 -fcV^ 



e '" " / (igQ-ik-{xi-Xj)^^-iki-ao _ ^+iki-do-^^^+ikj-ao _ ^-ikj-do-^ (4.27) 



27r JO 

where fcj = /cRi"^. Simplifying using equation |4.14|, we have 



W, 



lij — e 



[Jo{l\xi - Xj + aoi - aoj\) 



(4.28) 



+Jo{l\xi - Xj - aoi + aoj\) 
-Jo{l\xi - Xj - aoi - aoj\) 
-Jo{l\xi - Xj + + aoj\)] 

where aoi = Ri~^ao. 

• Case2: For just the diagonal window functions, equation |4.21| simplifies 



to 



Wku = {1/271) dct)\Bm{k)\\ 



(4.29) 
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The following are examples of equation [4.29| for simple beams and chop- 
ping strategies. 

Example 2.1. Single difference Gaussian beam. 
The diagonal window functions are: 

Wi = Bl{l)2[l - M2W,)] (4.30) 

where 

is the Fourier transform of an (unmodulated) symmetric Gaussian beam. 
Example 2.2. Double difference Gaussian beam. 
The diagonal window functions are: 

Wi = 52(/)[^ - 2Jo(/^c) + IM2Wc)] (4.32) 

Example 2.3. Single difference top hat beam. 
The diagonal window functions are: 

Wi = B^il)2[l - M2W,)] (4.33) 

where 

Bt{l) = [2MWt)]/iWt) (4.34) 

is the Fourier transform of a top hat beam of characteristic size 6t- 
Example 2.4. Double difference top hat beam. 
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The diagonal window functions are: 



Wi = B^S)[\ - 2^0 W + \U2iec)] (4.35) 



• Case 3: The diagonal window functions for a general modulation pat- 
tern, without assuming the flat sky appoximation are given by 

Wlmm' = B\l) MmsMm's'Pl{cOs{es - 9,,)) (4.36) 

s s' 

where 6 is the chopper angle as a function of sample s and Mm is the ampli- 
tude of modulation m. B'^{1) is the spherical transform of the (unmodulated) 



beam map. The beam must be symmetric. Some examples of equation |4.36 
are given below. 

Example 3.1. Single difference Gaussian beam. 

The modulation is given by: 

+1 at e = -Oc (s = 1) 

M = I (4.37) 

[ -1 at e = +ec {s = 2). 

and the spherical transform of a Gaussian beam is given by 

S(/) =e-^'('+^)"'. (4.38) 
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Thus, from equation [4.36|, the diagonal window functions are 



Wl = Bl{l) ELi Ekl MmsMm',,Pi{cOs{e, - ^,0) 

= Bl{i)[{+i){+i)Pi{cos{-e, + e,)) 
+{+i){~i)Pi{cos{-e,-e,)) 

+ (-l)(+l)PKcos(+^, + ^,)) (4.39) 

+{-i){-i)Pi{cos{+e,- e,))] 

= Bl{l)[Pi{cosm - Pi{cos{-2d,)) 
-Piicosi29,)) + PiicosiO))] 

Wl = B^M)2[1 - Pi{cos2e,)] (4.40) 



Example 3.2. Double difference Gaussian beam. 
The modulation is given by 



M = i 



-1/2 at e = -e^ (s = i) 

+ 1 at 9 = {s = 2) 
-1/2 at 9 = +9c {s = 3) 
and the diagonal window functions are given by 



(4.41) 



Wl = B^Mn - m{cos9,) + -Pi{cos29,)]. 



(4.42) 



Example 3.3. Four-point-chop Gaussian beam. 
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The modulation is given by 



M 



-1/4 


at 


9 = 




(. = 1) 


+3/4 


at 


e = 


-\0c 


(. = 2) 


-3/4 


at 


6 = 


+ \0c 


(. = 3) 


+ 1/4 


at 


e = 


+ |^e 


(. = 4) 



(4.43) 



and the diagonal window functions are given by 

-SIS 2 3 4 1 

Wi = Bl{l)[- - -Piicos-9,) + -Piicos-e,) - -Piicos2e,)]. (4.44) 

A summary of the diagonal window functions for simple chopping strate- 
gies is given in Table ||. The diagonal window functions in the fiat space 
approximation are similar to the ones calculated not assuming the fiat space 
approximation except that the Legendre polynomials Pi{cos6c) are replaced 



by Bessel functions Jq{Wc)- Figure ^ shows the diagonal window functions 
for single and double difference Gaussians with and without the fiat space 
approximation . 



modulation 



W, 



2- point 

3- point 

4- point 



2[1 - Pi{cos2ec)] 
[I ~ 2Pi{cose,) + lPi{cos2e,)] 
fPi{cosU,) + fP^cosl^c) - lPi{cos2e,)] 



Table 6 Summary of analytic expressions for diagonal window functions for 
simple chopping strategies. All are multiplied by B'^{1). The fiat space 
approximations are the same except the Legendre polynomials Pi{cos6c) are 
replaced by Bessel functions Jq{10c)- 
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Figure 30 Diagonal window functions with and without the flat sky approx- 

1 II 

imation. Shown in the top panel are 

W^'^ for single and double difference 
chops of Qc = 0.7° and a Gaussian beam of width a = 0.212°. Shown in the 
bottom panel are the residuals. 

4-2.5 Python V Window Functions 



For the PyV observing strategy, equation [4.21| becomes 



Wiij = (l/27r) J e-*^-(^'"^"^)5(fc)5*(fc)rf0 (4.45) 

where B{k) is the Fourier transform of the beam map for the given mod- 
ulation, Xi the position of field i, and k = /(cos 0, sin 0). These functions 



are computed for all pairs of fields and channels. Figure ^ shows the PyV 
diagonal window functions. A different modulated beam map is used at each 



declination because the total chopper throw, which is constant in azimuth 
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Multipole 1 

Figure 31 Diagonal window functions for PyV at a declination of —49°. The 
unapodized cosine modulation is plotted with a dashed line and the apodized 
cosine modulations are plotted with solid lines. 

degrees, subtends fewer real degrees on the sky at higher dechnations. 



4.3 KL Decomposition 



Karhunen-Loeve (KL) decomposition (see for example Tegmark et al. 
1996), also known as the signal-to-noise eigenmode method (Bond 1994), 
transforms the data into a basis in which modes of the data can be sorted 
by their signal-to-noise (S/N). 
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This technique is useful for testing the self-consistency of data sets, for 
examining modes with a particular S/N range, and for data compression. 
Calculation of the determinant in equation is an A^'^ process and cal- 



culation C^^d is an A^^ process (see for example the backslash command 
in MATLAB), where N is the number of data points. By eliminating the 
modes which do not contribute to parameter estimation, (i.e. the ones with 
low S/N), KL decomposition can compress data sets with a minimal loss of 
information, speeding up the data analysis process typically by 3 orders of 
magnitude. 

The steps for KL decomposition are as follows: 

1. Let CN' = + C^. 

2. Calculate (C^ )^^/^, for example by Cholesky decomposition. 

3. Calculate the matrix 

M = (CN')-i/2cT(c^')-*/2_ (4_4g) 

4. Calculate R, the matrix of eigenvectors of M. 

5. Sort R by eigenvalue. The eigenvalues are a measure of signal/noise. 

6. The compression matrix B' is the first A^keep rows of B = R(C^')~^/^, 
where A^keep is the number of KL modes to keep. 
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7. Then the theory matrix in the KL basis is 

= (B')C^(B')*, (4.47) 

the data in the KL basis is 

4l = {B')d, (4.48) 
and the noise matrix in the KL basis is the identity matrix I. 

4.4 Constraint Matrices For Offsets 

The subtraction of offsets can be accounted for as an additional term in 
the covariance matrix (e.g., Bond et al. 1998). For PyV, we must take into 
account the subtraction of the chopper synchronous offset and the ground 
shield synchronous offset. 

In both cases, the subtracted data in field i is given by: 

where /' is the mean of the offset vector which has been subtracted from 
field i. Both constraint matrices thus have the form 

CS = ^{Mj). (4.50) 
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The chopper synchronous offset is constant over a file and the same chop- 
per offset has been subtracted from fields i and j, thus /* = and the 
constraint matrix takes the form: 



(«:i/iVmcs) (4.51) 



for fields taken with the same set of files and zero for fields not taken with 
the same set of files. iVgies is the number of files which observed the fields. 
Even if the offset subtracted is not perfectly estimated, taking ki to be large 
ensures that we have no sensitivity to modes of the data that could have 
come from the chopper. An image of the chopper constraint matrix is shown 



in Figure 32 



In the case of the ground shield synchronous offset, /* ^ The ground 
shield constraint matrix is: 

Cf = K,{fJ,). (4.52) 

for fields i and j taken with the same set of files, and is zero for fields not 
taken with the same set of files. For each field fi is known from the data 
reduction process (section 3.5). There is a different ground shield constraint 
matrix for each modulation. Again, Hi2 is taken to be large. An image of the 



ground shield constraint matrix is shown in Figure 33 



We will refer to the total constraint matrix for PyV as C*^ = C*^^ + C*^^ 
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Figure 32 Image of chopper offset constraint matrix. C'"^{1, 1) is in the 
bottom left corner of the image. The matrix is 690 x 690 for the 345 fields 
and 2 channels in the reduced, co-added data. Points 1:345 correspond to 
channel 45 and points 346:690 correspond to channel 12. 
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Figure 33 Image of ground shield constraint matrix for modulation 2. 

is in the bottom left corner of the image. The matrix is 690 x 
690 for the 345 fields and 2 channels in the reduced, co-added data. Points 
1:345 correspond to channel 45 and points 346:690 correspond to channel 12. 



CHAPTER 5 
SINGLE MODULATION ANALYSIS 



The angular power spectrum is first estimated from tlie PyV data set by 
considering eacli modulation separately. An analysis which includes cross- 
modulation correlations and which estimates the angular power spectrum 
simultaneously in a series of /-space bands is discussed in Chapter 6. This 
single modulation analysis is simpler computationally and conceptually than 
the cross-modulation analysis. The reduced PyV data set contains 690 pixels 
X 8 modulations = 5520 data points. Computing the likelihood is consider- 
ably faster if 8 690x690 matrices rather than a 5520x5520 matrix is used. 
Additionally, the cross-modulation noise was difficult to model, as will be 
discussed in Chapter 6. 



5.1 Power Spectrum Constraints 

A fiat power spectrum is one for which C = (Z(/+l)Q/27r) is constant. For 
each of the eight modulations, we compute the likelihood (£) as a function 



of C = (/(/ + 1)Ci/2tt) following equation |4.2|. Limits are obtained from 



the likelihood curves following the prescription of Ganga et al. (1996), by 
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integrating L over C^/^ starting with the most hkely value of C^^^ and shcing 
at equal values of L until 68% of the total area under the curve is reached for 
Icr limits or 95% of the total area under the curve is reached for 2cr limits. 



An example likelihood curve and resulting limits are shown in Figure p| 




50 100 150 

(1(1+1)C/27t)'/^ (a^K) 



200 



Figure 34 Getting limits from a likelihood curve. The likelihood curve and 
Icr limits for modulation 6 are shown as a typical example. The shaded area 
is 68% of the total area under the curve. 



The CMB power spectrum is shown in Figure |3^ and is given in Table 
0. The band powers shown in Figure ^ are calculated for each modulation 
separately. The error bars include statistical uncertainties only and do not 
include uncertainties in the calibration or beam size. The calibration uncer- 
tainty allows all band powers to shift by the same amount (i.e. the calibra- 
tion errors are correlated). Given the beam size uncertainty of approximately 
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mode 


k 




1 
2 


74li 


231^ 
2611 


3 

4 
5 
6 
7 


108+1? 
14011? 

203li 

r)QQ+40 


3111 

281^ 
5411? 

96111 

9iti 


8 


26 





Table 7 Angular power spectrum constraints from the single modulation like- 
lihood analysis. Band powers are in iiK. 

0.015°, the band power for a given modulation can move roughly by a factor 
of exp(±/(0.425)(0.015)(7r/180)), only a 3% effect at / = 200. The effective / 
of each modulation is given by 



I{Wi 



(5.1) 



where 



(5.2) 



The I range of each modulation is determined by the half-maximum points 
of {Wiy/'^ and is indicated as horizontal error bars. 



5.2 Noise Model 



Our noise model assumes the covariance between fields taken with dif- 
ferent sets of files is neghgible because of the chopper offset removal and 
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Figure 35 PyV angular power spectrum constraints from the single modula- 
tion analysis. The detections have la error bars and the upper limit has 2a 
error bars. The unapodized cosine modulation is plotted with an open square 
and the apodized cosine modulations are plotted with diamonds. The error 
bars include statistical uncertainties only and do not include uncertainties in 
the calibration or beam size. The / range of each modulation is determined 
by the half-maximum points of {WiY^"^. Low / values correspond to large 
angular scales and high I values correspond to small angular scales. CMB 
power is clearly rising from low to high I up to the sensitivity cutoff of PyV. 

because of the long time (at least 10 hours) between measurements. An 
analysis comparing the noise covariance estimated from data which had not 
yet been co-added over all cycles and the noise covariance estimated from 
data which had been co-added over all cycles indicates that PyV noise is 
dominated by detector noise. However, the long term drifts due to the at- 
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mosphere, which add to the variance as well as induce small correlations 
between fields taken with the same set of files, could be important, especially 
for power spectrum estimation. 

The noise covariance between fields taken with the same set of files is first 
estimated by taking the usual covariance on the co-added data: 

NmesCS = E « - d'M - d^) (5.3) 

JVfiles - i n=l 



where 



]^ Affile 



d' = ^T. dl (5.4) 

-'^ files n=l 

d^ is the data from field i, d^ is the data from field j and A^aies is the number 
of files in the set containing fields i and j. However, since there were typically 
only 100 files taken for each field, the sample variance on the noise estimate 
is ~ 1/(100)^/^, or 10%, which will severely bias estimates of band power. 
To obtain a better estimate of the noise, we averaged the variances for each 
set of files and then scaled the off-diagonal elements of the covariance to the 
average variance in a given set based on a model derived from the entire PyV 
data set. 

In order to derive a noise model for the single modulation analysis we 
examined histograms of Affiles for each modulation as a function of stare 
lag from to +16 for channels 12 and 45 and from -16 to +16 for the cross- 



covariance between channels 12 and 45. Examples are shown in Figures 36 
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and We also examined plots of the mean of C^A^fiies for each modulation 
as a function of stare lag. If the off-diagonal correlations were due entirely 
to the chopper offset subtraction and white noise, the noise matrix would be 



of the form: 



Affiles = (Si 



(5.5) 



for fields i and j in the same observing set and noise variance a^. As il- 



lustrated in Figure 38 C A^fUcs is not flat as in equation |5~5|, but has some 



taper, indicating that long term drifts, perhaps due to atmosphere, could be 
important. The cross-channel covariances are small but non-zero. A sample 



noise covariance matrix is shown in Figure 39 
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Figure 36 Histogram of Cat Amnios for a stare lag of zero, modulation 5, channel 
45. 
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Figure 37 Histogram of CAriVfiies for a stare lag of two, modulation 5, channel 
45. 
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Figure 38 CjviVfiics vs. stare lag for modulation 5, channel 45 as an example. 
Cjv -^fiies is not flat, but is tapered for non-zero lags. 
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Figure 39 Noise covariance matrix, C^, for modulation 2. C^{1, 1) is in the 
bottom left. The matrix is 690 x 690 for the 345 fields and 2 channels in the 
reduced, co-added data. Points 1:345 correspond to channel 45 and points 
346:690 correspond to channel 12. 

5.3 Internal Consistency Checks 

Several self-consistency checks were performed on each modulation of the 
data set, using our best estimate of the noise covariance. We performed these 
same tests using the preliminary noise covariance, which suffers from sample 
variance, and found that the consistency tests failed. The tests indicate that 
the data set used is self consistent and that our best estimate of the noise 
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mode 


x' 


1 


615 


2 


660 


3 


603 


4 


600 


5 


667 


6 


710 


7 


640 


8 


609 



Table 8 at best fit band power. Expected dof = 690 - 31 x 2 = 628 
covariance is indeed a good model for the noise. 



5.3.1 Tests 

The = d^C-'^d, where cT is the data vector and C = + + 
is the total covariance, is consistent with its expected value (the number of 
degrees of freedom) in all modulations. Table ^ lists for each modulation. 



5.3.2 l3 Tests 

Since the two channels do not observe exactly the same points on the sky, 
a direct test between them could not be done. Instead, the probability 
enhancement factor, /3, (Knox et al. 1998) is computed between data from 
the two channels for each modulation. 
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The value of (3 is given by 
]^ln\C\ + ]^^C-^d- ^Zn|Cii | - ^^Cir' - \ln\C22\ - ^^,022-% (5.6) 

where Cn is the total covariance for the first channel, C22 is the total co- 
variance for the second channel, and C is the total covariance for all of the 
pixels in each modulation. The expected value of (3 is given by 



</5>o=2^n(^^j^J (5.7) 

and the variance of P is 

<{P- Pof >o= Tr(wi2W2i) (5.8) 

where 

W2i^Cji(Cri + Cri)-^ (5.9) 

and C21 is the theoretical cross-covariance between the two channels. 

The value of /3 falls within the expected range for all modulations when 
the best estimate noise matrix is used. This is not our strongest test be- 
cause the signal-to-noise is weaker when each channel and modulation are 
considered separately. 

5.3.3 KL Decomposition 



Finally, the data set was transformed into the KL (signal-to-noise eigen- 
mode) basis. In that basis, C~^/^d should be Gaussian distributed with 
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(7=1 and a total area equal to the number of degrees of freedom. His- 
tograms of C~^/^(i*are consistent with Gaussian distributions for all of the 
modulations for the best estimate noise matrices and are inconsistent with 



Gaussian distributions for the preliminary noise matrices. Figures ^ and 



41| show histograms for one of the modulations using the best estimate noise 



matrix and preliminary noise matrix respectively. 
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Figure 40 Distribution of C~-^^'^d in the KL basis for modulation 7 using our 
best estimate noise matrix. The large spike at C~^^'^d = is due to modes 
which have been set to zero by the chopper constraint matrix. The histogram 
follows the expected Gaussian distribution. 
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Figure 41 Distribution of C^^^'^d in the KL basis for modulation 7 using a 
preliminary noise matrix. The large spike at C~^^'^d = is due to modes 
which have been set to zero by the chopper constraint matrix. The number 
of counts in the central region are systematically lower than the number in 
the expected Gaussian distribution and there are many high a points. This 
preliminary noise matrix also failed the P tests. 



CHAPTER 6 
CROSS-MODULATION ANALYSIS 

An analysis including the cross-modulation correlations is presented in 
this chapter. This analysis is not simply a trivial expansion of the single 
modulation analysis presented in Chapter 5, but requires more sophisticated 
techniques for the estimation of the noise model and of the hkelihood. Fol- 
lowing the procedure of Wilson et al. (1999), the likelihood is simultaneously 
estimated in a series of 8 bands from a map which was constructed including 
all of the cross-modulation correlations. The results are consistent with those 
from the single modulation analysis. However, since this analysis includes 
the correlations between modulations, the constraints on the power spectrum 
can be used directly for cosmological parameter estimation. 

6.1 Noise Model and Consistency Checks 

In order to derive a noise model for the cross-modulation analysis, we 
expanded the procedure of section 5.2 for the cross- modulation noise cor- 
relations. The value of = ^(C'^ -I- C^)~^^, where C'^ includes noise 
correlations between all modulations, was computed for each observing set 
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Figure 42 x^/ dof for each set using an expansion of the noise model of section 
5.2. = d{C^ + C*^)^^(i* and typical dof = 13 stares x 7 modulations - 
7 constraints = 84. The x^/dof are unacceptable, so a better model of the 
cross- modulation noise is necessary. 



and channel. The results are shown in Figure 42 and indicate that a better 



model of the cross-modulation noise is necessary. 

In general, the elements of the cross-modulation noise matrix can be 



expressed as 



(6.1) 



where the components of C are given by 



(6.2) 



and Mm is the timestream modulation vector for modulation m, describes 
the noise in the timestream as a function of chopper sample s and is a 
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model of the noise between fields i and j. Thus models for C'^ and are 
needed. For clarity, is a 128 x 128 matrix, C'^ is a 8 x 8 matrix, and 
is a 690 x 690 matrix for PyV. 

In the single modulation analysis we effectively assumed was the iden- 
tity matrix. For white noise with the chopper offset removed, C'^ would be 
of the form: 

C^mm' OC (5., - -^^M'^Mm' (6.3) 
'stares 

However, a single a does not characterize all modulations. Rather, a de- 
creases with increasing modulation, so a more sophisticated model for is 
necessary. 

Assuming only depends on chopper sample separation As = s — s', 
we can compute from the following function: 

f{As) = J2dJs+As. (6.4) 

s 

For example, the Cfj component of is given by /(As = 1). In order to 
compute /(As), a chopper synchronous offset is first subtracted from the raw 
data. Then /(As) is calculated for each channel, cycle and stare in a file. 



/(As) is then averaged over cycle, stare, and file. Figure |43| shows /(As) for 
each channel in one of the sets. 



is calculated for each set and channel following equation |6.2|. C 



matrices for channels 1 and 2 right and left are averaged, as are matrices 
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Figure 43 /(As) for set hb9. /(As) breaks around As = 6, as would be 
expected from the beam size. 

for channels 4 and 5 right and left, yielding matrices for channels 12 
and 45 in each set. As an example, for one set and channel is shown in 
Figure H^. 



In order to get a simple form for C^, we investigated the effect of ignoring 
the cross-channel correlations and the taper as a function of stare lag on the 



single modulation band powers. The results are shown in Figure 45. Ignoring 



the correlations between channels 12 and 45 and assuming the correlation 
between fields i and j come only from the chopper offset subtraction does 
not change the angular power spectrum significantly, so we assume is of 
that form for the cross-modulation analysis. 



Finally, since this noise model is derived from sample-sample, it is larger 
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Figure 44 for set hb9, channel 45. Elements that are more than 2 mod- 
ulations apart are relatively uncorrelated. 

than the corresponding noise derived from the co-added data by a factor 
of ~ 10^, so Cfj^^i must be normalized to the variance in the co-added 
data for each set. Since C'^ accounts for the relative normalization of all 
of the modulations, Cf^^^t must be normalized to the variance in only one 
modulation of the co- added data for each set. We normalize to modulation 8 



because it is least potentially affected by the ground shield. Figure E6| shows 



the / dof for each set using the final cross-modulation noise model. If 



I J mm' 

2 



is normalized to modulation 2 instead of modulation 8, the x depends on 
elevation (Figure H^. 



As another check on the noise matrix used in the cross-modulation anal- 
ysis, single modulation band powers were computed using the CjJ^^ compo- 
nents of C^j^j^,. They are consistent with the band powers given in Chapter 
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Figure 45 Angular power spectra obtained using 4 different noise models. Ig- 
noring the correlations between channels 12 and 45 and the taper of CNNfUes 
as a function of stare lag does not change the angular power spectrum sig- 
nificantly. 



5. 



6.2 Map Making 

The power spectrum is simultaneously estimated in several /-space bands, 
this time including all of the cross-modulation covariances, from a map con- 
structed from the modulated data. The power spectrum could have been 
estimated using covariance matrices of size 5520 x 5520 in each of the /-space 
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Figure 46 x /dof for each set using the final noise matrix, x = d{C + 
C^)~^c/* and typical dof = 13 stares x 8 modulations - 8 constraints = 96. 



The vertical scale of the plot covers a smaller range than that of Figure 42. 



bands, but transformation into a map basis first reduces the computational 
load. This technique was used in an analysis of the MSAM-I experiment 
(Wilson et al. 1999), for which the power spectrum estimated using the 
map is consistent with the power spectrum estimated directly from the full 
covariance matrices. 

In order to construct the map, first recognize that the data can be ex- 



pressed as 



d = BT + n 



(6.5) 



with noise covariance matrix < nn >= N. (The actual N used is N = 
C'^ + C^.) To obtain the underlying temperature field T, the matrix B 
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Figure 47 x^/dof for each set normalized to modulation 2. = d{C^ + 
C^)~^(i* and typical dof = 13 stares x 8 modulations - 8 constraints = 96. 
The declination of the sets roughly decreases with increasing set number. 
The vertical scale of the plot covers a larger range than that of Figure HB. 



must be inverted. This inversion is carried out by constructing the estimator 
T which minimizes the x^'- 



X = id-BT)N-^{d-BT). 



(6.6) 



We find 



(6.7) 



This estimator will be distributed around the true temperature due to noise, 



where N, the noise covariance matrix for the map, is given by 



N=<{f- f){f - f) >= (b'^N^^b) \ 



(6.8) 
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Table 9 Angular power spectrum constraints from the cross-modulation anal- 
ysis. Band powers are in fiK. 

6.2.1 Power Spectrum Constraints 



This map can be analyzed in the same manner as the modulated data, 
with the advantage that the signal covariance matrix is much simpler to 



compute than using equations 4.3 and 4.21. In the map basis, the signal 



covariance matrix simplifies to 



2/ -I- 1 



(6.9) 



Indeed, one way to think of a map is that it is the linear combination of the 
data for which the signal (and therefore its covariance) is independent of the 
experiment. The noise covariance for the map (equation |6.8|) accounts for all 
of the experimental processing. 

The results of the likelihood analysis are shown in Figure |8 and are given 
in Table |9|. 



82 



150 



100 



+ 50 



50 100 150 200 

Multipole 1 



250 



300 



Figure 48 Angular power spectrum constraints from the cross-modulation 
analysis. Flat band power, (/(/ + 1)Ci/2t[Y^'^ is simultaneously estimated 
in 8 Z-space bands including cross-modulation theory and noise covariances. 
The error bars include statistical uncertainties only. 

The angular power spectra obtained using the cross-modulation and sin- 



gle modulation analyses are consistent, as illustrated in Figure At first 
glance, the band power in the fourth modulation of the single modulation 
analysis is higher than that of the fourth band of the cross-modulation anal- 
ysis. However, the Knox minimum variance filter (Knox 1999) shows that 
the fourth modulation of the single modulation analysis has significant sen- 
sitivity to lower / scales. The Knox filter is a more complete representation 
of angular sensitivity than just the diagonal window functions, because it 
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Figure 49 Comparison of angular power spectra from single and cross- 
modulation analyses. Points obtained from the single modulation analysis 
are denoted by gray squares and points obtained from the cross-modulation 
are denoted by black diamonds. 



contains information from the full (off-diagonal) covariance matrix. 



6.2.2 CMB Images 

The map can also be Weiner-filtered to produce a realistic image of the 
sky. A Weiner-filtered map of the main PyV region and the S/A^ in the map 
are shown in Figure ISO. The unfiltered map, not the Weiner-filtered map, 



was used for power spectrum estimation. 
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6.3 Comparisons With Other Experiments 

6. 3. 1 Power Spectrum Comparison 

Flat band powers calculated from the subset of PyV data in the Pylll 
region of sky are consistent with the Pylll flat band powers to within the 
uncertainties. Given the calibration and statistical uncertainties, the lowest 
/ PyV modulation agrees with the the smallest scale COBE measurements. 



6. 3. 2 Comparing Py V and Pylll Maps 

Using the technique of section 6.2.1, maps of the sky are reconstructed 
from the Pylll data. The Pylll map is compared to the PyV map in the 



same region of the sky in Figure Since Pylll has higher S/N than PyV, 
features which appear in the PyV map should also appear in the Pylll map, 
but features which appear in the Pylll map may not necessarily appear in the 
PyV map. In order to compare modes of the experiments which are similar 
to each other, both maps are transformed into the KL basis and only modes 
with S/N > 2 are used. It is not necessary to Weiner filter the maps (which 
removes noise) if only the highest S/N modes are used. Structures found in 



the PyV map of Figure are also evident in the Pylll map of Figure El 



implying that Pylll and PyV are consistent with each other. 
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Figure 50 Weiner-filtered CMB and S/N maps for the main PyV region. The 
unfiltered map was used for power spectrum estimation. 
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Figure 51 A comparison of PyV and Pylll maps. Structures found in the 
PyV map are also evident in the Pylll map, implying that Pylll and PyV 
are consistent with each other. 



CHAPTER 7 
FOREGROUNDS 

The PyV data are cross-correlated with several foreground templates in 
order to set limits on possible foreground contamination. The templates used 
are the Schlegel et al. (1998) 100 micron dust map, which is based on IRAS 
and DIRBE maps, the Haslam et al. (1974) 408 MHz survey (synchrotron), 
and the PMN survey (point sources). Each foreground template map is 
smoothed to PyV resolution, pixelized and modulated according to the PyV 
observation scheme. 

Two templates are created for the PMN survey. We call one PMN, which 
is converted to STcmb using the spectral indices given in the survey. The 
other we call PMNO, which is converted to a flux at 40 GHz assuming a 
flat spectrum extrapolated from the flux measurement at 4.85 GHz. The 
assumption of a flat spectrum is conservative in that it is likely to over- 
estimate the flux at 40 GHz. Neither case is correct, since spectral indices 
have not been measured for all of the sources, in which case a flat spectrum is 
assumed, but we do know that some of them are not flat, so a flat spectrum 
will be inappropriate. The two cases cover a reasonable range of possibilities. 

For each modulation and foreground template, a correlation coefficient 
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and uncertainty are calculated following de Oliveira-Costa et al. (1997): 

a = i^;al = ^^ (7.1) 

—* — * 

where df and dp are the modulated foreground and PyV data respectively 
and C = + + (Figure 

A weighted mean and uncertainty over all of the modulations for each 



foreground are given in Table [T^. In all cases there is no clear detection of 



foregrounds. The RMS of each modulation of each foreground was calculated 



and then multiplied by the corresponding la error bar in Table 10 in order to 



estimate an upper limit on the foreground contribution to CMB band power. 



The limits on contributions from foregrounds are given in Table |Tl| and are 
at least ~ lOx smaller than the measured CMB bandpowers. 

As a test on both the map-making procedure and the foreground estima- 
tion procedure, a map of the dust was constructed from the dust data which 
had been modulated and pixelized according to PyV observing strategy. The 
original dust map was recovered. 

If the diffuse morphology of the sky is not constant as a function of wave- 
length, then these templates do not reveal all of the foreground contamination 
and more could be hidden in the PyV data. A combined analysis of Pylll at 
90 GHz and PyV at 40 GHz would constrain the foreground contamination 
further. 
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Figure 52 Correlation coefficients and uncertainties for each foreground and 
modulation. 
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Table 10 Correlation coefficients and uncertainties for weighted means. 
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Table 11 Upper limits on foreground contribution. All units are ^K. 
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Figure 53 Upper limits on foreground contribution. All units are ^K. 



CHAPTER 8 
CONCLUSIONS 

The goal of the Python V season was to increase the sky coverage and 
the range of observed angular scales of the Python experiment thereby con- 
straining the CMB angular power spectrum. 

The PyV experiment fully samples 598 deg^ of the microwave sky and 
constrains the CMB angular power spectrum in the angular scale range 40 ^ 
/ ^ 260. The measurements pass internal consistency checks, show little 
contamination from foreground radiation, and are consistent with previous 
Python and COBE results. 

The results of experiments previous to PyV collectively imply a rise in 
the angular power spectrum, but the PyV experiment definitively shows a 
rise in the spectrum from from large (/ ~ 40) to small (/ ~ 200) angular 
scales. Figure ^ summarizes the constraints on the angular power spectrum 
by experiments previous to PyV and highlights the PyV results. 

The PyV data are inconsistent with a fiat spectrum, and instead show 
a rising spectrum, implying the acoustic scenario of structure formation de- 
scribed in the introduction is correct. The PyV data give confidence in results 
of estimation of cosmological parameters, such as Q, from all previous CMB 
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Figure 54 CMB anisotropy measurements as of Feb 1999. The points are the 
same as those in Figure |^. The PyV measurements are overplotted. 

experiments because the rise in the spectrum could no longer be a calibra- 
tion artifact. Since reionization of the universe would damp power at smaller 
angular scales (Hu and White 1995), the PyV data can rule out models of 
homogeneous reionization with Zreion ^ 75. Earlier work on cosmic defect 
models suggests there would be no peak in the CMB angular power spec- 
trum due to defects (Allen et al. 1997, Pen et al. 1998), in which case the 
PyV data would disfavor defect models. However, there is still disagreement 
in the predictions, with later work (Albrecht et al. 1999, Wu et al. 1998) sug- 
gesting there could be a peak due to defects. The PyV spectrum is shown in 
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Figure 55 PyV angular power spectrum and COBE-normalized standard 
CDM model (density Q = 1.0, Hubble constant Hq = 50 km/s/Mpc, baryon 
density fib = 0.05, and cosmological constant A = 0.0). Only statistical 
uncertainties are included in the PyV points. 

reference to a COBE-normalized standard cold dark matter (sCDM) model 



in Figure |5^. The PyV spectrum has lower power at large angular scales 
and a steeper slope than typical sCDM models (x^/dof ~ 2.5 for the model 
shown). Models with a higher baryon density and/or a cosmological constant 
have a steeper slope than the model shown. While future CMB experiments 
will estimate cosmological parameters to high precision in the coming decade, 
the robust rise in the Python V spectrum already tightly constricts the range 
of viable models of structure formation. 
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